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1. Introduction 

A defect-free crystal is a macroscopic quantum system with spatial periodicity. In 
principle, only nonlocal observables are probed by plane waves, eg. with diffraction or 
spectroscopy techniques. This is attested by the marked localization in reciprocal space 
of Bragg peaks as well as discrete phonon states. However, since the seminal work of 
Slater [TJ, the nonlocal nature of protons is commonly overlooked in crystals containing 
bistable O-H- ■ -O hydrogen bonds, when the coexistence of two configurations, say 
"L" for 01 — H • • ■ 02 and 11 R" for 01 ■ ■ ■ H — 02, is conceived of as classical disorder. 
This leads to a conflict of interpretation consisting in that a statistical distribution of 
protons should destroy the translation periodicity of the crystal, whereas infrared and 
Raman spectra evidence extended states at the Brillouin-zone center (BZC) obeying the 
symmetry-related selection rules of the crystal. This conflict has been tolerated because 
the nonlocal nature of protons cannot be evidenced by x-ray diffraction. Even with 
neutrons, a Bragg peak analysis is not conclusive, for statics and dynamics cannot be 
unravelled effectively from the unit-cell structure factor [2]. 

The same conflict of interpretation arises for crystals containing hydrogen bonded 
centrosymmetric dimers, such as carboxylic acids (3j H], or potassium hydrogen 
carbonate (or bicarbonate, KHCO3) [SI El El IH]- The thermal interconversion 
between tautomeric configurations, LL and RR, is supposed to be correlated with 
a reorganization of the single and double CO bonds, along a reaction path across a 
multidimensional potential. This is commonly rationalized with rigid proton-pairs in 
an asymmetric double-well coupled to an incoherent thermal bath mimicking a liquid- 
like environment. It is supposed that "phonon-induced incoherent tunnelling" should 
prevail at low temperatures, whereas stochastic jumps should take the lead at elevated 
temperatures. By contrast, vibrational spectra are consistent with nonlocal observables 
represented by normal coordinates. There is no evidence of any symmetry breaking due 
to local disorder and the density-of-states remains discrete at any temperature [9| [TO], [TTj . 

Nowadays, the local or nonlocal behaviour of bistable hydrogen bonds can 
be distinguished at advanced pulsed neutron sources, with a time-of-flight Laue 
diffractometer measuring a much broader range of reciprocal space than any reactor- 
based four-circle instrument. The scattering function measured at large momentum 
transfer values (Q = kj — ky, where kj and kj are the initial and final wave vectors, 
respectively) contains information on the scale of spacial correlations between L and R 
sites (~ 0.5 A). Moreover, such an instrument is unique to observing diffuse scattering, 
because signals in addition to Bragg peaks are effectively measured across the whole 
range of reciprocal space for each neutron pulse. It is thus possible to determine whether 
or not coherent scattering is affected by the coexistence of L and R configurations at 
various temperatures. 

This has been exemplified with the crystal of KHCO3 made of hydrogen bonded 
dimers, (HCOj) 2 [321 EH El EH]. Below 150 K, only those proton sites (L), 
corresponding to a unique configuration for dimers (LL), are occupied. The neutron 
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scattering function reveals enhanced Bragg peaks exclusively at those particular de- 
values corresponding to the reciprocal sublattice of protons. Furthermore, the enhanced 
peaks are superposed to a prominent continuum of incoherent scattering, cigar-like 
shaped along Q- vectors perpendicular to dimer planes. By contrast, there is no evidence 
of similar enhanced features for the isomorphous derivative KDCO3. To the best 
of our knowledge, the dramatic increase of both coherent and incoherent scattering 
intensities for the sublattice of protons, not observed for deuterons, cannot be explained 
with any conventional model based upon Slater's localized particles. Alternatively, a 
quantum theory leads to different macroscopic quantum correlations (entanglement) for 
indistinguishable protons (fermions) or deuterons (bosons). The enhanced features can 
be thus explained by the spin symmetry of the "super-rigid" fermion states for which 
phonons are forbidden. By contrast, deuteron states are represented by phonons with 
no spin-symmetry. 

Above 150 K, the less favored sites along the hydrogen bonds (R), at A x rj 0.6 
A from the L sites, are progressively populated, while the space group symmetry 
and the unit cell parameters remain unchanged. At 300 K, the population ratio is 
~ 18 (RR) : 82 (LL). The enhanced features observed continuously from cryogenic to 
room temperatures show that the quantum correlations are temperature independent. 
The coexistence of LL and RR configurations is consistent with a quantum superposition 
of proton states perfectly insulated from decoherence by the other degrees of freedom 
p~3j HI] . Then, scattering functions, vibrational spectra, and quasi-elastic neutron 
scattering (QENS) data, can be reinterpreted consistently with nonlocal observables and 
quantum interference arising from the superposition of macroscopic tunnelling states. 

Structural transitions of second order [P1\ja < — > C2/m) occur at T c h ~ 318 K 
(KHCO3) or T c d Pd 353 K (KDCO3). Above T c , dimers split into two centrosymmetric 
entities, symmetric with respect to the plane (a, c). It is quite remarkable that the 
enhanced features of the scattering function, as well as proton dynamics and quantum 
interference, do not show any discontinuity at T c [TJ]. There is no evidence of 
classical disorder above T c . It transpires that quantum correlations arise from the space 
periodicity of centrosymmetric dimers, irrespective of the space group symmetry. 

Here, we examine whether the enhanced features survive, or not, in a mixed 
isotope crystal of KH p Di_ p C03. In this case, static or dynamical disorder are clearly 
differentiated and the existence of distinct crystallographic domains of KHCO3 or 
KDCO3 is excluded. A static distribution of mutually exclusive particles should destroy 
the quantum entanglement and the enhanced features should vanish. Alternatively, if 
isotopes merge into separable macroscopic states, such that quantum correlations across 
the sublattice of protons survive, the enhanced features should be observed, although 
with lower contrast compared to KHCO3. To the best of our knowledge, such a nonlocal 
approach has never been examined in depth, in this context. 

In a preliminary report [IE] , the structure of KH0.76D0.24CO3 at 300 K was refined 
from Bragg peaks measured with the four-circle diffractometer 5C4 at the Orphee 
reactor of the Laboratoire Leon Brillouin [17J. There was no evidence of any symmetry 
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breaking: the space group [P2\/a) is identical to that of KHCO3 or KDCO3 at 
the same temperature (Fig. [1]). A structural model consistent with the proposed 
quantum theory was represented by four sublattices, B.ll, Hrr, Dll, D_rr, with weights 
PhQllh 1 PhQrrh 1 PdQlld ■> PdQrrd > respectively: 



The bound coherent scattering lengths are bu = —3.741 fm or bp = 6.671 fm. The 
refined positional and thermal parameters (R ~ 3%) were identical for protons and 
deuterons whose sites are indiscernible, to within experimental precision and thermal 
ellipsoids. The estimated concentrations were pu = 0.76, pr> = 0.24, and the relative 
weights were such as S = PhQllh ~\~ PhQrrh ~i~ PdQlld ~\~PdQrrd = 0.99(4). There was no 
evidence of any significant decrease of Bragg peak intensities into Laue monotonic diffuse 
scattering due to a statistical distribution of isotopes. Alternatively, for a sublattice of 
paired up nonseparable protons and deuterons, the averaged scattering length should 
be: 



The model is more artificial because p is no longer a free parameter and it is less 
satisfactory (R « 4%). However, it cannot be definitely excluded. Consequently, the 
question as to whether the sublattices of protons and deuterons are separable entities 
remains largely open. 

The conclusions of the present work are based upon the neutron scattering function 
of the mixed isotope crystal presented in Sec. El The enhanced features show that the 
quantum correlations of the pure crystal are not destroyed by partial deuteration and 
Raman spectra are consistent with separable proton and deuteron states. In Sec. El the 
quantum theory [HI [15] applied to the mixed crystal accounts for macroscopic proton 
and deuteron states, each of which being a superposition of LL and RR states. A 
detailed interpretation of the enhanced features is proposed. In Sec. HI we show that 
alternative explanations are inappropriate. 

2. Experimental 

2.1. The neutron scattering function 

The mixed isotope crystal (~ 30 mm 3 ) measured in the previous diffraction work [T6] 
was cut from a twinning-free parent crystal (~ 20 x 15 x 3 mm 3 ) obtained from a water 
solution (^ 75% H 2 0, 25% D 2 0). This parent crystal was loaded in the vacuum chamber 
of the time-of-flight Laue diffractometer SXD, at the ISIS pulsed neutron source [181I19] . 
The temperature was (300 ± 1) K. We measured 8 orientations, each for « 8 hours, with 
{b) || (&*) either parallel or perpendicular to the equatorial plane of the detector bank. 
The reduced data were concatenated in order to compute the scattering function (Fig. 
[2]). Previous studies of crystals with similar sizes have shown that multiple-scattering 
effects are not important. However, Bragg peaks show severe extinction. A standard 



Ph[qllhHll + orrhRrr] + Pd[qlld^ll + ^rrd^rr]- 



(1) 



b = 0.76b H + 0.24b D ss -1.242 fm, 



(2) 
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analysis is consistent with the structure and isotope contents previously determined, 
but the modest precision (R « 10%) shows that measurements conducted at different 
sources are complementary. 

The interpretation of the scattering function is based upon the crystal structure 
made of centrosymmetric dimers lying in (103) planes with O- • -O bonds parallel to each 
other. Proton/deuteron coordinates are x,y,z, (Fig. [Q, Q x ,Q y ,Q z , are projections of 
Q, U xx , Uyy, U zz , are the dominant temperature factors. Dotted lines in Fig. [1] enhance 
the network of double lines of proton/deuteron sites, either parallel to (6), separated by 
A x « 0.6 A, or parallel to x, with A y 2.2 A. The spatial periodicity of these grating- 
like structures is D x « a/ cos42° « 20.39 A or D y = b. The distance between (x,y) 
planes is D z w 3.28 A. A previous analysis of rigid body motions has shown that there 
is no short-range correlation between protons and carbonate entities p2] and inelastic 
neutron scattering (INS) spectra show that proton dynamics are separated from heavy 
atoms and virtually dispersion-free [201 EI] • 

In Fig. [21 cuts of the scattering function parallel to (a*, c*) show spots of intensity, 
corresponding to Bragg peaks, and, at large Q-values, enhanced streaks along Q z 
perpendicular to the trace of dimer planes parallel to (Q x ,Qy), with Q y \\(b*). These 
streaks are not an artefact of the instrument because: (i) the detectors were calibrated 
with vanadium; (ii) the streaks are measured by different detectors for each of the 8 
orientations of the crystal; (iii) spurious signals eventually due to an anomalous detector 
should appear along different trajectories. 

The maps in Fig. [2] are similar to those of KHCO3 at the same temperature |13j . 
apart from the weaker intensity of diffuse scattering due to the smaller concentration 
of protons. On the other hand, they are markedly different from those of KDCO3 
for which enhanced features are totally extinguished [12]. The enhanced diffraction 
peaks correspond to nodes of the reciprocal sublattice of protons/deuterons and there 
is no obvious alternative sublattice with similar parameters. They occur when k = 
Qy/b* = ±2Tm y /{b* A y ) fh ±n y x 2.57. For k = 0, in-phase scattering by double 
lines of sites parallel to x gives peaks at Q x = ±2ir/A x = ±(10.00 ± 0.25) A -1 and 
Q z m n z x 1.92 A -1 . At k = ±3, n y m 1, anti-phase scattering along Q y and Q x give 
peaks at Q x = ±(5.00 ± 0.25) A^ 1 , Q z rts n z x 1.92 A -1 . For k = ±5 (n y ns 2) and 
k = ±8 (n y ~ 3), the patterns of enhanced peaks are similar to those at k = and 
k = ±3, respectively. 

A cut along the streak at Q x = (10.00 ± 0.25) A" 1 in the plane k = (Fig. 0J) 
shows enhanced Bragg peaks on the top of a continuum along Q z . Because there is no 
evidence of disorder along z, this continuum is essentially due to incoherent scattering 
[12]. It compares favorably with the Gaussian profile I w exp — Q 2 Z U ZZ , with U zz ~ 0.04 
A 2 [16J. Compared to the intensity J at Q x = 0, Q z = 0, the ratio Io/ho ~ 2, 
with substantial uncertainty, is much greater than expected from conventional models 
according to which incoherent scattering by protons should be proportional to the Debye- 
Waller factor DWF = exp — (Q X U XX + Q^Uyy + Q 2 Z U ZZ ), while incoherent scattering by 
other nuclei should be negligible. With U xx ~ 0.035 A 2 [16], the expected intensity ratio 
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should be Iq/I[ q ~ 30. The observed intensity (/ 10 ) suggests that U xx is dramatically 
reduced for these Q x , Q y - values. 

Fig. H] (top) shows a cut at k = 0, parallel to a*, through the reciprocal plane 
at —(5.0 ± 0.2) A -1 off-a*. The enhanced Bragg peak observed at —(8.5 ± 0.2) A -1 
(spotted by an arrow) hides the underlying incoherent scattering signal. Alternatively, 
a cut at —(4.2 ± 0.2) A -1 off-a* is free of Bragg peaks (Fig. HI bottom). It reveals 
the sharp profile of incoherent scattering centered at —(9.0 ± 0.2) A" 1 (spotted by an 
arrow). Compared to the enhanced Bragg peak, the shift in position by ~ —0.5 A -1 is 
consistent with the slope of the Q z axis. The half-width at half-maximum, A a * 0.4 
A -1 , gives the half-width along Q x , Aq x m A a */ cos 42° « 0.3 A" 1 , similar to that of the 
Bragg peak (^ 0.2 A -1 ). Consequently, U xx vanishes exclusively if Q x matches a node of 
the reciprocal sublattice of protons (see below). Otherwise, incoherent scattering is not 
observed, presumably because it is dramatically depressed as exp — (Q 2 X U XX ). Similarly, 
cuts around k = (not shown) suggest that U yy vanishes sharply outside the plane. 

2.2. Separation of proton and deuteron states: Raman scattering 

In Fig. [51 the Raman profile of the stretching mode z/OD for the pure crystal is red- 
shifted by ~ 600 cm -1 with respect to the z/OH profile of KHCO3 and the spectrum of the 
isotope mixture is a mere superposition of these profiles. There is no visible frequency 
shift for each profile and, therefore, there is no significant dynamical correlation between 
isotopes. By contrast, if protons and deuterons were paired up, a single z/0(H p Di_ p ) 
profile should appear at an intermediate frequency, shifting smoothly from z/OH (p = 1) 
to z/OD {p = 0). This is clearly excluded. The spectra accord with separable proton 
and deuteron extended states consistent with the mixture of sublattices ([T|). 

In addition, the O- • -O bond lengths are significantly different for KHCO3 (Rh = 
2.58 A) or KDCO3 (R D = 2.61 A) and #0.76 = 2.59 A for KH0.76D0.24CO3 is seemingly 
a weighted average [16]. Isomorphous sublattices of protons and deuterons exclude 
a statistical distribution of O- • -O distances. Otherwise, the z/OH frequency should 
vary sig nificantly with p, for Az/OH/Ai^ « 5000 cm^A- 1 [22J. The absence of 
frequency shift, or band broadening, suggests that a mixed crystal can be represented by 
isomorphous sublattices pfKHCOs] and (1— p)[KDC03], whose structures and dynamics 
are independent of p. These sublattices are not resolved by diffraction. 

3. Theory 

Admittedly the enhanced features are an evidence of quantum correlations, but different 
mechanisms have been envisaged: symmetry related entanglement [121 HSl [EH 115] : 
quantum exchange [23]; spin-spin interaction [23]. They will be compared. 

The Raman spectra suggest that the Hamiltonian for a mixed crystal can be 
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partitioned as 



Hp 



— p(Hh + %kco 3H + Ch,KC03h) 

+ (1 — p)v~Lg + %kco 3D + Cd,kco 3D ) + ^ 



(3) 



"Hh, "Hd, represent the sublattices of H + , D + ; "Hkco 3H > %kco 3D , those of heavy nuclei; 
Ch,kco 3 h> ^d,kco 3D ! couple light and heavy nuclei. The sublattices are bound through the 
electronic term Hep- There is no dynamical coupling between the KHCO3 and KDCO3 
subsystems which are supposed to be separable: any state vector of 7i p is a tensor 
product of state vectors of the subsystems. Furthermore, the adiabatic separation of 
Hn and "Hkco 3H > or and "Hkco 3D , leads to factorable state vectors for the subsystem 
of bare protons (fermions), or deuterons (bosons) [23j 1211 125] . The separation of proton 
and deuteron states is, therefore, consistent with the "super selection rule" for fermions 
and bosons [26J . 

The mixed crystal is described with jV 0! Nb, N c (A/" = N a NbN c ) unit cells labelled 
j, k, I, along crystal axes (a), (6), (c), respectively. The sublattices are described with the 
same unit cell containing 4p KHC0 3 and 4(1— p) KDC0 3 entities. Proton/deuteron sites 
are organized in two dimer entities indexed as j, k, I and j', k, I, respectively, with j = j' . 
For each dimer, dynamics are represented by normal coordinates, either antisymmetric 
{aw} or symmetric {a sr } (r = 1,2,3), which are linear combinations of displacements 
along x,y,z, respectively. The coordinates of the centers of symmetry, {a 0r }, {a:' 0r }, 
are (x ,y ,z ), (x' Q , y' , z' ). For planar dimers, z = z' = 0. Configurations LL and 
RR correspond to opposite signs for x and x' . Although proton and deuteron sites 
may be not strictly identical, their dynamics are formally represented by the same set of 
normal coordinates. The wave functions in the ground state are represented by gaussian 
functions ^Qj kl (a ar ) and ^ s jkic( a sr — V^aorc), with C = li LL n or "RR". 

A key issue of the theory is to clearly differentiate sites and states. Sites 
refer to local positional parameters, whereas each single state is extended over all 
indistinguishable sites. Conversely, the probability density (occupancy) at any site 
comprises contributions from every single state. The nonlocal observables (normal 
coordinates and conjugated momentums) are not bound to any particular unit cell. 
We should emphasize that nonlocality due to the lattice periodicity is not related to 
any physical interaction. It is energy free. This is totally different in nature from the 
long-range quantum exchange in a supersolid, like the crystal of He [27]. Obviousy, 
the proposal by Keen and Lovesey [23] that nonlocality could arise from non vanishing 
exchange integrals is incorrect, even for the nearest sites separated by ~ 2.2 A [28J. 

3.1. Macroscopic states 

It is customary to represent lattice dynamics with plane waves (phonons). For deuterons 
in the ground state, this can be written as: 

, N c N b N a 




1=1 k=l j=\ 
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where k is the wave vector; Ljjy = ja + kh + Ic is a lattice vector; a, b, c, are the unit 
cell vectors. The wave functions for a unit cell, E 0D j k[i; c = ®oDjkic + ^oDj'kic with q = 
"+" or "— ", is a linear combination of the wave functions for centrosymmetric sites: 
®0Djkic({a r}) = Yl^0Djki( a ar)%Djkic( a sr-V2a 0rC ) and likewise for Q Dj'kic{W 0r }). 

r 

For centrosymmetric dimers, there is no permanent dipole, so long-range correlations 
through inter-dimer interaction are negligible. Phonons are virtually dispersion-free 
[201 EI]. 

|D£.)c is a single-deuteron state with an occupancy number (4AT) -1 per site. This 
state cannot be factored with eigen states for local coordinates (xijki, x 2jki ■ ■ •)• It i s 
entangled in position. (By definition, a pure state of an arbitrary number of particles 
is said to be entangled when, because of its indivisible nature, it cannot be described 
by treating each particle separately.) Likewise, the state in momentum space is also 
entangled. In addition, this state meets the symmetry postulate for bosons stipulating 
that the sign should not change when any two dimer sites, say j, k, I, and j", k", I", or 
any two sites within dimers, are interchanged [29J. 

Alternatively, eigen states for fermions should be antisymmetric. An interchange 
of lattice sites j, k, I, and j", k", I", should transform a plane wave into the anti-phase 
wave and superposition of these waves should annihilate the amplitude. Hence, phonons 
are forbidden, except those matching the lattice periodicity, for which any interchange 
is neutral. Such states represent in-phase dynamics in each unit cell, namely ^ = "+" , 
and in-phase dynamics for all unit cells, namely k ■ Lh = modulo 2tt, where Lh is 
any vector of the sublattice of protons. This macroscopic quantum effect excluding 
elastic distortions has been named "super-rigidity" [13]. In addition, eigen states must 
be antisymmetrized with respect to interchange of proton sites within dimers, what 
is equivalent to changing the signs of the coordinates for the center of symmetry, say 
{a 0r } — > —{a 0r }. To this end, the wave function of a dimer, supposed to be non planar 
for the ease of notation, is rewritten as a linear combination of the wave functions for 
permuted sites as: <3>' 0HjklCT = ^ ll2 [Qmjkic{{o^r}) + rQ Hjkic(-{ocQr})} (and likewise 
for j'), with t = "+" or "— ". Then, antisymmetrization requires a spin-symmetry, either 
singlet-like (5* for r = +) or triplet-like (T for r = — ). Finally, with the translation 
invariant wave function of the unit cell, ^oHjki+Cr = ®oHjkiCT Jr ®oHj'kiCTi protons states 
can be written as: 

N c N b N a 

\ R X)c = -7j7 2222z2\ E oHjki+c+)\ s )^ 



1 N c N b N a 
\ E -)c = -rr7^2^2^2\ E 'oHjkl+C-)\ T )- 



(5) 



N 1=1 k=l j=l 

For each configuration C, the ground states are entangled in position, momentum, and 
spin. The singlet-like and triplet-like states are degenerate because the spin symmetry 
is energy-free. The proposal by Sugimoto et al. [24J that the spin symmetry could be 
due to spin-spin interaction is incorrect because the actual coupling between nearest 
spin sites (~ 10 4 Hz ~ 10 -6 K) is far too weak to give rise to long-lived correlations. 
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For an isolated crystal, symmetry-related quantum-correlations are decoherence- 
free. Thanks to the adiabatic separation, there is no phonon-induced decoherence and, 
therefore, no transition to the classical regime, before melting or decomposition. In 
practice, some states can be transitorily disentangled by external perturbations, such 
as thermal radiations, but this is counterbalanced by spontaneous re-entanglement 
via decay to the ground state. In a standard gaseous atmosphere, the density of 
the surroundings is ~ 1000-fold smaller than that of the crystal, so the impact 
of environment-induced decoherence to measurements is invisible. A source of 
disentanglement is the thermal population of excited proton states (for example 7 a 
OH and 7 S OH at « 1000 cm -1 ), as the a — s splitting destroys the spin-symmetry, but 
this is marginal. 

Because the spin symmetry holds for both LL and RR ground states, for which 
there is no a — s splitting, the sublattice at 300 K can be represented by a mixture 
comprising a coherent superposition of super- rigid states for protons (E]), on the one 
hand, and, on the other, a coherent superposition of phonon states for deuterons (BJ: 

\h + (t)) = Uq'LLH(T)\K)LL + gWCniH+W]; 

M(T)) = hq' LLD (T)\D^) LL + q' RRD (T)\Dl) RR ]; (6) 
with |<Zi£#(T)| 2 = qLLniT), etc, in accordance with (TjQ). 



3.2. Neutron scattering 

The partial differential nuclear cross-section can be written as [30] 

d^[£ = pE^i^/i^i^wi^i^)! 2 ^ + 3» - ( 7 ) 

1 hf 

The summation runs over initial states weight pi, and final states |/), Ef. The 

energy transfer is fkv and dQ is an element of solid angle. The interaction potential 
between the incident neutron and the target sample, V'(r), comprises the interaction 
potential at every site. For diffraction, \i) = \f), = Ef, hu> = 0. Neutrons are 
diffracted by each single state. Because of the adiabatic separation, neutrons diffracted 
by proton or deuteron states are free of any contribution from heavy atoms. Moreover, 
neutrons diffracted by proton states probe exclusively the interaction potential at the 
proton sites and likewise for deuteron states, irrespective of whether or not these sites 
are identical. There is no Laue monotonic diffuse scattering. In general, V"(r) does 
not have the periodicity of the lattice, for the potential will depend upon the nuclear 
position and spin orientation. However, this is not the case for the super-rigid states 
which match exactly the periodicity of the crystal. 

Upon momentum transfer to the target proton states (jBJ), the final state is, in 
general, such that kf-Ln ^ modulo 2tt, so elastic distortions are generated. They 
are represented by ellipsoids at H/D sites in Fig. [TJ The super-rigid states are probed 
without induced distortions when 

Q • L H = modulo 2tt. (8) 
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Then, neutrons probe the spin-symmetry and the coherent cross-section is the total 
nuclear cross-section, oh ~ 82.0 b, that is ~ 45 times greater than the coherent cross- 
section for regular Bragg peaks, o~h c ~ 1-8 b, or « 15 times greater than a^c ~ 5.59 
b, or 3-fold greater than o"kco 3 c ~ 27.7 b for heavy atoms. These scattering events 
emerge from the diffraction pattern with spectacular contrast, even at large |Q|-values 
as they are not attenuated by the Debye- Waller factor. Besides, they are temperature 
independent. 

For the P2\/a structure, z is not parallel to an axis of the sublattice of protons (Fig. 
HJ, so the matching condition (JS} is never realized along Q z . The super-rigidity along z is 
necessarily destroyed by measurements and intensities are proportional to exp —U ZZ Q 2 Z . 
On the other hand, since Zq = z' = for planar dimers, the spin-symmetry holds only in 
the (x, y) planes, so the differential cross-section for enhanced diffraction can be written 
as [Hj 



Tj Tf 



K N b N c 

Yj J2 Y ex ViQzlD z {[expiQy (kD y - A y /2) 
j=i k=i i=i 

+ TiTf exp iQ y (kD y + Ay/ 2)] ^ 
x [expiQ x (jD x - A x /2) + TiT f expiQ x (jD x + A x /2)}} 2 
x exp(-U zz Q 2 z ) x 6{u), 

where {• • -} 2 takes into account indistinguishable scattering events, either in-phase 
(\S) — > \S) and |T) — > \T)), or anti-phase (\S) — > \T) and |T) — > \S)). Enhanced 
diffraction by (x, y) planes is anticipated at Qu y = ±27m y /A y , Qhx = ±2im x / A x for 
n y even (in-phase scattering) and Qhx = ±2^(277,^ + 1)/A X , for n y odd (anti-phase 
scattering). Then, enhanced Bragg peaks occur at Q z = 2im z /D z , with intensities 
proportional to exp —U ZZ Q 2 Z . 

Otherwise, if Q z ^ 2nn z /D z , the super-rigid (x,y) lattice is probed without 
distortion for the same (Qhx, QijjJ-values, but the spin-symmetry is destroyed. 
Incoherent scattering is observed and the intensity is proportional to 

a m exp -(Q 2 Z U ZZ )5(Q X - Q Hx + Q y - Qh v ), (10) 

hence the cigar-like shape. Compared to a nonrigid lattice, the enhancement is merely 
due to U xx = U yy = 0. The observed intensity ratio Iq/I\q ~ 2 is at variance with the 
expected ratio of ~ 1.5, but this discrepancy could be representative of errors in the 
background correction. 

Figures [20 can be thus rationalized with three different kinds of scattering events: 
(i) regular Bragg diffraction by the crystal lattice; (ii) enhanced diffraction by the fully 
entangled sublattice of protons; (iii) incoherent scattering by the super-rigid planes of 
protons. 
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4. Discussion 

The proposed theory is based upon fundamental laws of quantum mechanics applied 
to the crystal in question: the structure is periodic; dimers are centrosymmetric; 
indistinguishable protons are fermions; indistinguishable deuterons are bosons. It leads 
to macroscopically entangled states and, in the special case of protons, to super-rigidity 
and spin-symmetry with observable consequences. This theory is consistent with a large 
set of experimental data (neutron diffraction, QENS, INS, infrared and Raman) and, to 
the best of our knowledge, there is no conflict with any observation. There is, therefore, 
every reason to conclude that the crystal is a macroscopic quantum system for which 
only nonlocal observables are relevant. 

Because quantum entanglement emerges naturally from the crystal structure, there 
is no reason to introduce any further interaction. From the viewpoint of economy 
(Ockham's razor), alternative theories based upon quantum exchange [23] or spin-spin 
interaction |24J are less satisfactory. Moreover, these correlations are very unlikely to be 
significant and it is doubtful that these short-range interactions, essentially confined to 
dimers, could account for macroscopic correlations, especially at elevated temperatures. 

Protons are unique to demonstrating quantum entanglement, because they are 
fermions and because the very large incoherent cross-section can merge into the total 
coherent cross-section. No other nucleus can manifest such an increase of its coherent 
cross-section. The enhanced features can be, therefore, unambiguously assigned to 
protons, in accordance with their positions in reciprocal space. They are evidences 
of macroscopic quantum correlations which have no counterpart in classical physics. 

According to the theory, the pair correlation function for protons represented by 
ellipsoids in Fig. [TJ are due to measurement induced phonons. The structure free 
of induced distortions should comprise Dirac functions for protons and ellipsoids for 
deuterons. The spatial wave functions for protons (jSJ) describe the distribution of the 
center of mass of the proton states with respect to that for all other atoms and the 
uncertainty principle holds at each site of the perfectly correlated sublattice. 

It is widely believed that macroscopic objects are, in principle, governed by 
the laws of quantum mechanics and classical physics is an emergent phenomenon at 
an appropriate limit when the system interacts with an environment that might as 
well consist of its internal degrees of freedom [31] . The crystal of bicarbonate is 
a counterexample that is never amenable to classical physics, because (i) quantum 
entanglement is as robust as the crystal; (ii) the adiabatic separation insulates the 
proton states from every decoherence mechanism; (iii) the energy gap for excited proton 
states is much greater than kT. 

A different approach places the emergence of classicality on the impossibility 
of distinguishing individual quantum levels with coarse-grained measurements [32J. 
Transposed to diffraction, the interaction potential in ((Tj) is replaced by the mean 
potential, V(r) = ^2Pi(i\V(r)\i), that is then Fourier transformed into the unit-cell 

i 

structure factor [30]. Eigen states and nonlocal observables are thus translated into 
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classical variables for particles located at definite sites and neutrons scattered by various 
nuclei can interfere, irrespective of whether or not these nuclei are involved in the same 
eigen states. Because of the opposite signs of the scattering lengths for protons and 
deuterons, interferences could build up in such a way that some Bragg peaks would 
be more intense for KHCO3 than for KDCO3, or vice versa. However, these enhanced 
peaks are not expected to single out the reciprocal sublattice of protons and the structure 
factor does not account for the enhanced incoherent scattering along Q z in any way. In 
addition, the rigid-body analysis does not provide any evidence of correlation between 
protons and heavy atoms [12J. There is, therefore, experimental support to conclude 
that the assumption of classicality emerging from coarse-grained measurements is only 
a convenient approximation hiding the quantum nature of the crystal. Nevertheless, 
quantum effects embedded in the Bragg peaks transpire as the mean scattering length 
((2]) is less satisfactory than the mixture of distinct sublattices (JTJ , although the difference 
is not large enough to be firmly conclusive. 

5. Conclusion 

Neutron diffraction and Raman spectra show that the mixed isotope crystal is not 
a statistical mixture of KHCO3 and KDCO3 entities. The crystal is a macroscopic 
quantum system made of separable matter waves over the isomorphous sublattices of 
KHCO3 and KDCO3. The quantum theory shows that deuteron (boson) states can be 
represented by phonons whereas antisymmetrized proton (fermion) states demonstrate 
super-rigidity and spin-symmetry. This macroscopic quantum entanglement is energy 
free and decoherence free. It is as robust as the crystal structure and the emergence 
of the classical regime at elevated temperatures is forbidden. The main outcome of 
the present work is that macroscopic quantum effects for protons are not destroyed 
by deuteron states. Despite their different masses, isotopes are not defects leading to 
localization, possibly because they are isoelectronic. 

The quantum theory of diffraction by proton states allows us to distinguish 
different events. First, in the most general case, momentum transfer destroys the 
quantum entanglement and regular Bragg peaks are observed. Second, for those 
particular Q-values matching the nodes of the reciprocal sublattice, the super-rigidity 
and the spin-symmetry are probed without any measurement induced symmetry 
breaking. The incoherent cross-section merges into the total coherent cross-section 
and the corresponding peaks show markedly enhanced intensities featuring macroscopic 
quantum entanglement. Because of the anomalous Debye- Waller factor, these peaks are 
observed at large Q-values and they are temperature independent. Third, for those Q- 
values matching only the nodes of the reciprocal planes of protons, the spin-symmetry is 
destroyed and cigar-like shaped rods of incoherent scattering by the super-rigid planes 
are observed in between the enhanced peaks. Clearly, the unit-cell structure factor is 
not appropriate to account for these enhanced features. 

This work presents one single case of macroscopically entangled states on the scale 
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of Avogadro's constant in a mixed isotope crystal at room temperature. The quantum 
theory suggests that such macroscopic quantum effects should be of significance for 
many hydrogen bonded crystals. Hopefully, unforeseeable consequences of the quantum 
view advocated in this work will appear in the future. 
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Figure 1. (Color on line) Schematic view of the crystalline structure of 
KH0.76D0.24CO3 at 300 K. Top: projection onto (a, b) planes. Bottom: projection 
onto (a, c) planes. The solid lines represent the unit cell. Proton and deuteron sites 
are indiscernible within thermal ellipsoids. The dashed lines emphasize the sublattice 
of protons/deuterons. The arrows point to the sites populated at low temperature. 
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Figure 2. Diffraction patterns of KH0.76D0.24CO3 at 300 K in (a*,c*) planes. The 
arrows point to the ridges of intensity for momentum transfer parallel to (z) and 
perpendicular to the trace of dimer planes (dashed lines), which contain the (x) 
direction, as defined in Fig. [TJ The (y) direction parallel to (&*) is perpendicular 
to (a*,c*) planes, k = Q v /b* . 
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Figure 3. Cut of the scattering function along Q z in the plane k = 0, at Q x = 
(10± 0.25) A -1 , after background correction. The Gaussian profile exp(— Q 2 z u 2 z ), with 
u z — 0.04 A 2 is representative of the Debye- Waller factor for protons determined by 
neutron diffraction in Ref. [16j . 




Figure 4. Cuts of the scattering function along a* in the plane k = 0, at — (5.0± 0.2) 
A _1 (top) and —(4.5 ± 0.2) A^ 1 (bottom) off-a*, after background correction. The 
arrows point to the enhanced Bragg peak (top) and to the diffuse scattering (bottom) . 
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Figure 5. Raman spectra of powdered crystal of KHCO3 (top), KH0.5D0.5CO3 
(middle), and KDCO3 (bottom) at 300 K. The dashed lines are guides for the eyes. 
Spectra were recorded with a triple-monochromator DILOR XY® equipped with an 
Ar+ laser emitting at 4880 A. The spectral resolution was w 2 cm -1 . The samples 
were sealed in glass tubes. 



